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ABSTRACT 


This paper presents computational fluid dynamic studies of the inlet turbulence and wall curvature effects on the flow 
steadiness at near wall surface locations in boundary layer flows. The time-stepping RANS numerical solver of the 
NASA Glenn-HT RAN S code and a one-equation turbulence model, with a uniform inlet turbulence modeling level of 
the order of 10 % of molecular viscosity, were used to perform the numerical computations. The approach was first 
calibrated for its predictabilities of friction factor, velocity, and temperature at near surface locations within a transitional 
boundary layer over concave wall. The approach was then used to predict the velocity and friction factor variations in a 
boundary layer recovering from concave curvature. As time iteration proceeded in the computations, the computed 
friction factors converged to their values from existing experiments. The computed friction factors, velocity, and static 
temperatures at near wall surface locations oscillated periodically in terms of time iteration steps and physical locations 
along the span-wise direction. At the upstream stations, the relationship among the normal and tangential velocities 
showed vortices effects on the velocity variations. Coherent vortices effect on the velocity components broke down at 
downstream stations. The computations also predicted the vortices effects on the velocity variations within a boundary 
layer flow developed along a concave wall surface with a downstream recovery flat wall surface. It was concluded that 
the computational approach might have the potential to analyze the flow steadiness in a turbine blade flow. 


INTRODUCTION 

The Gas Turbine Branch of NASA Glenn Research 
Center has been developing the NASA Glenn-HT CFD 
Code for the prediction of gas turbine blade surface heat 
transfer. The Code used time-stepping RANS 
computational methods and two equation turbulence 
models to perform the surface heat transfer computations. 
Wall curvature and inlet turbulence play important roles 
in the flow development over a turbine airfoil. This paper 
presents the approach and results of using the RANS 
numerical solver and a One-Equation Turbulence Model 
to model the effects of inlet turbulence level and wall 
curvature on flow steadiness at locations near wall 
surface. 

The author previously proposed 1 to use the RANS solver 2 
and a version of One-Equation Turbulence Model 3 to 
analyze the vortices effects on flow properties with 
surface coolant injection into a boundary layer flow. The 
author then used the approach to perform a preliminary 
study 4 of vortices effects on flow with the influence of 
concave wall curvature. The previous studies 1,4 showed 
that the computational approach, as time iteration 
proceeded, could predict the experimental friction factors 
and span-wise velocity variations in boundary layer flows 
with vortex structures. The author reports here a study of 
using an existing computational approach 1 ' 4 to capture the 


inlet turbulence and wall curvature effects on the flow 
steadiness at near wall surface locations in boundary layer 
flows. The approach was first calibrated according to the 
experiments 5 of transitional boundary layer on a heated 
concave wall. 



Fig 1 Vortices induced inflection velocity profiles 5 
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The experiments measured the velocity and temperature 
variations to investigate the Taylor-Gortler vortices 
occurrence within the boundary layer flow. When vortices 
occur, they produce cross-span and cross-stream 
inflection point velocity profile, Figure 1, and the flow 
can be steady, transitional, or unsteady. The surface 
temperature measurements 5 , Figure 2, showed vortices 
also induced cross-span temperature variations. This 
author used some of the existing experimental results to 
calibrate the numerical method. 



Fig 2 Vortices induced span-wise temperature profile 5 


As a preliminary step to test the approach for analyzing 
the steadiness of turbine blade flow, the approach was 
used to predict the velocity and skin friction variations 
from other experiments 6 . The experiments measured 
velocities and friction factors of a boundary layer over a 
concave wall surface followed with a flat recovery wall 
surface. The velocity measurements showed vortices 
establishment and vortices effect on boundary layer flow 
development. This author tested the computational 
approach to predict some of the velocity measurements. 

Time stepping numerical computations were performed 
with specified inlet turbulence parameter, downstream 
pressure, computational domains and computational grids. 
The computed friction factors and the velocities and 
temperatures at near surface locations were compared 
with their values from the existing experiments. The three 
dimensional velocity and temperature variations were 
investigated in terms of time iteration steps as the 
computed friction factors converged to their experimental 
values. The study concluded that this computational 
approach could be improved for unsteady turbine blade 
boundary layer flow analysis. 


NOMENCLATURE 


Dimensionless Symbols 


Cf 

friction factor 

n 

number of time iteration steps 

N 

distance from the wall surface 

Pd 

exit static pressure 

S 

stream-wise coordinate 

T 

static temperature 

v n 

normal velocity along N direction 

V s 

tangential velocity along S direction 

z 

span-wise coordinate 

X 

velocity profile wave length 

Xi 

inlet turbulence parameter 

Reference Conditions 

h 

0.375 ft 

Pr 

14.7 psia 

Tr 

535° R 

v r 

58 ft/sec 

Pr 

density at T r and P r 

V 

molecular kinematic viscosity at T r 

p 

molecular viscosity 


RANS COMPUTATIONS 


The author used the RANS numerical solver of the NASA 
Glenn-HT code to perform the numerical computations. 
The solver used a time-stepping numerical scheme 7 to 
solve the three-dimensional, compressible, turbulent time- 
dependent Navier-Stokes equations. An explicit four stage 
Runge-Kutta numerical scheme 2 and fourth order cell- 
center spatial discretization scheme were used for the 
present computations. The Numerical solutions of the 
Navier-Stokes equations and the One-Equation turbulence 
model were advanced in time using a local time step 2 , At. 
This time step accounted for both the convective ( At c ) 
and diffusive contribution (AG) according to At = Co(At c x 
At d )/(At c + At d ). The mathematical expressions for At c , At d 
and stability control could be found in previous work 2 . Co 
was a constant and was similar to Courant-Friedrichs- 
Lewy number. This author assigned a constant value of 1 
to C 0 in the computations. The computational domains 
and corresponding computational grid configurations, 
turbulence model, and boundary conditions were 
described in the following. 

Computational Domains and Grids 

The physical domain used for the method calibration was 
chosen according to the test section in the experiments 5 . 
The test section was a rectangular turning duct. The 
height of the test section was 0.375 ft (h in this paper) and 
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test section had a high aspect ratio (width/height = 5.6). 
The radius of the concave wall surface was 3.18 ft (8.5h). 
The duct lower wall was adjusted to develop a boundary 
layer flow along the concave wall surface. The length of 
the concave wall surface extended over an angle of 79 
degrees and the wall length was 11.75h. The 
measurements were made within a small center region of 
the upper concave wall. Thus, this author chose a physical 
domain, Figure 3, corresponding to the upper half of the 
turning duct for the computations. 


>« 



Fig 3 Schematic of domain, 0-3-4-5-6-7-2-1, 
for method calibration. 


Although the domain size along the Z direction was small 
(0.175h), the size enclosed the domain of velocity 
measurements in the experiments 5 . Gridgen software was 
used to generate the computational grid. A single block 
structured computational grid, with grid point dimensions 
of 173, 73, and 57 along S, N and Z directions, was first 
used in the computations. The first grid point along the N 
direction was at a distance of 0.00005 from the wall 
surface. This distance was kept the same along the S 
direction. The grid point spacing increased gradually to 
0.03 at the domain edge (N = 0.5). The grid point spacing 
along S direction was stretched from 0.059 at the inlet (S 
= 0) to 0.18 at the exit (S = 11.75). The grid point spacing 
was uniform along Z direction. To study the grid effect on 
the computations, the computations were repeated with a 
grid point dimension of 205, 61, and 57. The first grid 
point along N direction was at N = 0.0001 and the grid 
point spacing was stretched to 0.03 at the N = 0.5 
location. The grid point spacing along the S direction was 
stretched from 0.059 at the inlet to 0.09 at the exit. The 
grid point spacing was uniform along Z direction. 


0.35h 



Fig 4 Schematic of domain, 0-3-4-5-6-7-2-1, 
for method application. 


For the method application, the physical domain, Figure 
4, was chosen according to the geometry of the test 
section in existing experiments 6 . Different from the 
domain in Figure 3, the concave wall surface was 
followed with a flat wall surface. The flat wall surface 
was tangent to the concave wall surface at the wall 
junction. The length of the flat wall was 12.74h. The 
domain size along the Z direction had been increased to 
0.35h to enclose the domain of velocity measurement 6 . 
The computational grid point dimensions were 405, 61, 
and 37 along the S, N, and Z directions. The first grid 
point, along the N direction at all S stations, was 0.0001 
from the wall surface. The grid point spacing was 
stretched to 0.03 at the domain edge (N = 0.5). Along the 
S direction, the grid point spacing was 0.03 at the inlet 
and was stretched to 0.09 at the exit. The grid point 
spacing was uniform along the Z direction. 
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Turbulence Model 

A one-equation turbulence transport equation, derived 
from an original turbulence model 3 , was used to model 
the eddy viscosity. The capability of this model to 
account for the flow transition with vorticity at a tripping 
point prompted the use of the model. The turbulence 
transport equation was written in a form 1 suitable for its 
implementation in the numerical scheme used in the 
Glenn-HT RANS code. The mathematical expressions of 
the turbulence model were summarized in the 
APPENDIX section of the paper. It was further assumed 
that each surface grid point along the span-wise direction 
beneath a field point acted as an individual tripping point. 
The transition term in the turbulence model was the 
summation of the transition terms due to vortices at all 
span-wise surface grid points. Reynolds analogy was used 
to relate the turbulence contribution in the energy 
equation. 

Boundary Conditions 

The surfaces of the computational domains (Figures 3 and 
4) consisted of inlet (0-1-2-3) and exit (4-5-6-7) planes, 
wall surface (0-3-4-5), and symmetry planes (2-3-4-7, 0- 
1-6-5 and 1 -2-7-6). Since a fourth-order cell-center 
scheme was used for the spatial discretization in the 
numerical schemes. The numerical computations required 
flow properties at two phantom cell centers immediately 
outside of the domain surfaces. These properties were 
specified with the following approaches. 

It was assumed that the cross flow velocity components 
were zero at the inlet phantom cell centers. Their static 
pressures and mass flow rates were set to the first interior 
cell centers’ static pressures and mass flow rates 
computed from previous time iteration. The density and 
total energy were then computed from the pressure and 
velocity with isentropic flow relationships. A uniform 
turbulent eddy viscosity parameter y, (APPENDIX), of 
the order of 0.1, was imposed at the inlet phantom cell 
centers. According to equations (2), (3), and (4) in the 
APPENDIX, this range of y; value gave a uniform inlet 
eddy viscosity level which was much smaller than the 
molecular viscosity. 

A zero normal pressure gradient was assumed at the wall 
surface and the linearization approach 8 was used to 
calculate the density values at the phantom cell centers 
next to the wall surface. Non-slip conditions were used to 
set the velocity and turbulent eddy viscosity at the 
phantom cell centers. The constant wall temperature 
condition was used in the computations and the 
temperature was 4 % larger than the room temperature 5 ' 6 . 


Flow symmetry was assumed at the domain surfaces (2-3- 
4-7 and 0-1 -6-5) along the Z direction and at the domain 
surface (1 -2-7-6) at N = 0.5. The flow properties at the 
phantom cell centers next to these surfaces were set to 
their values at the interior cell centers adjacent to the 
boundary surfaces. 

A uniform static pressure was imposed at the exit plane. 
The other mean flow properties at the down stream 
phantom cell centers were extrapolated 8 from their 
interior values and the exit static pressure. The present 
author used a downstream static pressure level, p d = 0.96, 
for the model calibration and a downstream static 
pressure, p d = 0.93, for the model application. It was 
assumed that the turbulent eddy viscosity values at the 
downstream phantom cell centers were the same as eddy 
viscosity value at the first interior cell center next to the 
exit plane. 

RESULTS AND DISCUSSIONS 

The computational procedures were performed for four 
different cases. The domain, grid point dimensions, inlet 
turbulence parameter, and downstream pressure for each 
case were listed in Table 1. 


case 

grid points 
S, N,Z 

domain 

turbulence 

X 

Pd 

1 

173,73,57 

fig 3 

0.050 

0.96 

2 

173,73,57 

fig 3 

0.075 

0.96 

3 

205,61,57 

fig 3 

0.075 

0.96 

4 

405,61,37 

fig 4 

0.125 

0.93 


Table 1 Parameters for computational cases 


The method calibration was based on the accuracy to 
predict friction factors from experiments 5 . Three 
computational cases (1,2, and 3) were performed to study 
the effects of grid point dimensions and inlet turbulence 
parameter, y; on the friction factor computations. With 
experience from the method calibration, a computational 
grid and an inlet turbulence parameter, y; = 0.125 were 
chosen for the method application (case 4). As time 
iterations advanced to the order of 10,000, the computed 
friction factor converged to its experimental value. The 
computed friction factor and the normal and tangential 
velocities and static temperature variations at locations 
near the wall surface were compared with their 
measurements to access the accuracy of the computational 
methods. As friction factors converged to their 
experimental values, the velocity and temperature 
variations in terms of the time iteration step and span- 
wise locations were studied for the steadiness of boundary 
layer flows. 
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Method Calibration 


It was assumed that the tangential velocity changed 
linearly from the wall surface to the first cell center 
location. The surface friction was computed with |iV s /N 
corresponding to the first cell center flow conditions. Skin 
friction was non-dimensionalized with reference 5 
dynamic pressure 0.5p r V, 2 to get friction factor, C f . C f at 
Z ~ 0(near the symmetric plane) from the computations 
were plotted in figure 5. 



S 

Fig 5 Friction factors at symmetric plane, z ~ 0 

Experimental results 5 of Cf at four stream- wise stations (S 
= 3.12, 5.35, 7.68, and 9.91) were also plotted in this 
figure to compare them with the computed C f at those S 
stations. For case 1 (Xi = 0.05) and case 2 ( Xi = 0.075), 
the computed Cf values, with a grid point dimensions of 
173, 73, and 57, seemed to converge at time iteration 
steps > 9000. The computations, with Xi = 0.05, predicted 
high Cf at S < 3.5 stations and very small Cf at S > 3.5 
stations. The computations, with Xi = 0.075, calculated 
higher Cf values than the case with Xi = 0.05, at S > 3.5 
stations. Cf at S = 3.12, 5.35, 7.68, and 9.91 from the 
above two cases did not agree well with the experimental 
values. With grid dimensions of 205, 61, 57 and Xi = 
0.075 (case 3), computed Cf values seemed to converge at 
n > 13700 and the computed Cf values at S = 3.12, 5.35, 
7.68 and 9.91 agreed very well with their experimental 
values. However, computed C f values were unsteady and 
oscillated when the time iteration steps were greater than 
13700. Figure 6 showed the C f variations (case 3) along 
the S direction at Z ~ 0 and three time iteration steps, n = 
13700, 14200, and 14650. These results showed that Cf at 
4.0< S < 6.5 and n = 14650 retained their values at n = 
13700 after significant reductions at n = 14200 time 
iteration. 



s 


Fig 6 Friction factors at z ~ 0 for 3 iteration steps 



Z 

Fig 7 Span-wise variations of V s at S = 3.12 and 5.35 


Existing experiments 5 measured the velocity at N = 0.002 
or 0.003 and Z < 0.1 locations at S = 3.12 and S = 5.35 
stations. Computed V s values (case 3) corresponding to 
these locations and n = 14650 were plotted in Figure 7 
and they were compared with their measured values. Both 
computations and measurements showed similar span- 
wise velocity variations. At S = 3.12 station, the V s from 
the computations and experiment had approximately the 
same wave length, X. However, the computations 
predicted a small amplitude at S = 3.12. The computations 
showed an increase in the amplitude of the V s variation at 
S = 5.35 but the experiments showed small amplitude of 
V s variation at S = 5.35. 
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Fig 8 Span-wise variations of V n and V s at S - 3.12 



z 



Fig 9 Span-wise variations of V„ and V s at S = 5.35 
and 7.68. 


Some of the computational results (from case 3) of 
tangential and normal velocity components at the first cell 
centers were studied for the patterns of the velocity 
variations over the concave wall surface. The Z direction 
variations of V s and V n , at two time iteration steps (n = 
14200 and 14650) and S = 3.12 station were shown in 
Figure 8. At Z ~ 0.05 and 0.13, V n was high and directed 
towards the wall surface at n = 14650 and this V n was 
associated with higher V s than its value at n = 14200. At Z 
= 0.088, V n was directed away from the wall surface at 
these two time steps. V n was high at n = 14650 and was 
associated with low V s . These relationships among the 
velocity components agreed with the effect of upwash and 
downwash on the velocity variations (Figure 1) at the near 
surface locations The Z direction variations of V s and V n 
at S = 5.35 and n = 13700, 14200 and 14650 and at S = 
7.68 and n = 14650 were shown in Figure 9a and b. V s at 
S = 5.35 oscillated at these three time iterations (Figure 
9a) and V s at n = 14650 retained approximately its value 
at n = 13700 with a different value at n = 14200. The 
effect of the V n direction (upwash and downwash) on the 
V s value appeared at n = 14650 time step. At S = 7.68 
station, V n (Figure 9b) changed its direction periodically 
at locations along the Z direction. However, the V s 
variation along the Z direction (Figure 9a) was not well 
correlated with the change in the V n direction (Figure 9b). 



z 

Fig 10 Span-wise Temperature variations at S = 5.35 
And S = 7.68 

The span-wise variations of the static temperatures, 
obtained from the computations, at the first cell centers at 
S = 5.35 for n = 13700, 14200, and 14650 and at the first 
cell centers at S = 7.68 for n =14650 were plotted in 
Figure 10. The temperature at S = 5.35 with n = 1370 and 
14650 and at S = 7.68 with n = 14650 all showed 
periodically span-wise variations. The temperature 
variation agreed well with patterns (Figure 2) of wall 
surface temperature obtained from the liquid crystal 
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technique 5 . The temperature at S = 5.35 also oscillated in 
terms of time steps at n = 13700, 14200, and 14650. 

Thus, the computations (case 3), with Xi = 0.075 and grid 
point dimensions of 205, 61, 57, predicted experimental 
observations of the near surface velocity and temperature 
within a boundary layer developed along a concave wall 
surface. The concave wall surface induced vortices at the 
upstream stations (S < 3). The calculated vortices induced 
upwash and downwash at near surface locations. 
Downwash increased the tangential velocity and upwash 
decreased the tangential velocity. The vortex broke down 
at the mid stream-wise stations (S = 5.35). The 
association of vortices and the velocity components 
disappeared at downstream stations and unsteady three 
dimensional variations in the velocity and temperature 
occurred at the near wall surface locations. 


Method Application 

As a preliminary step to analyze the turbine blade flow 
with the method, it was used to compute (case 4) the near 
wall surface velocity and friction factors of flow over a 
concave wall with a downstream recovery flat wall 
surface 6 . 



Fig 11 Stream-wise variations of C f at symmetric 
plane, z ~ 0 


Computed Cf at Z ~ 0 (near symmetric plane) and all S 
stations were shown in Figure 11. Cfat five different time 
iteration steps within 9400 < n < 9900 were presented in 
this figure. Cf values at 6 stream-wise stations from the 
experiment 6 were also plotted in this figure to verify the 
accuracy of the Cf calculations. The computed Cf at n = 
9400, 9630, and 9900 agreed well with the experimental 
C f at S = 7.5, 9.9, and 12.2. The computed C f values at the 
flat wall surface (S > 12.5) were lower than their 


experimental values. However, the computed Cf at the flat 
wall surface was about 2.5xl0" 3 . This value was close to 
the friction factor of a fully developed turbulent boundary 
layer flow along a flat wall surface. At all S stations, the 
computed Cf value was almost the same at n = 9400, 
9630, and 9900. At n = 9510 and 9750, C f values at the 
concave wall surface were higher than their values at n = 
9400, 9630, and 9900. This indicated that V s at the near 
wall surface locations oscillated in terms of the time 
iteration step. As an example, the computed (case 4) V s 
and V n at S = 3.35, Z = 0.077, 0.175, and 0.272 within 
9300 < n < 10000 were plotted in Figure 12. Both V s and 
V n oscillated periodically in term of time iteration steps. 
V n also changed its direction periodically and showed 
downwash and upwash effects on V s . 



n 

Fig 12 V n and V s variations for 9300 < n < 10000 


The computed V„ and V s values at N = 0.00005 locations 
along the Z direction at S = 3.35 and n = 9490, 9575, 
9663, and 9775 were shown in Figure 13(a) and (b). At n 
= 9490 and 9663, V n directed away (upwash) from the 
wall surface and had approximately the same values at n 
= 9490 and 9663. At n = 9575 and 9775, V n directed 
towards (downwash) the wall surface and V n values were 
very close at all Z locations. V s values at 0.08 < Z < 0.27 
oscillated at these four time steps. V n and V s at 0.08 < Z < 
0.27 showed that positive V n was associated with low V s 
and negative V n was associated high V s . This V n and V s 
relationship showed the vortices effects on the span-wise 
velocity variation. 
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z 

Fig 13 V„ and V s correlations for 9490 < n < 9775 


Computed V s values at n = 9600, N = 0.0032, 0 < Z < 
0.35 and S = 5.5 and 7.7 were plotted in Figure 14 and 
they were compared with the velocity measurements from 
the experiment 6 . The computed V s values at S = 5.5 
agreed well with the velocity measurements at Z < 0.15 
locations. The computed Vs oscillated along the Z 
direction, especially at S = 7.7 station, and the 
measurements showed only slight oscillation. However, 
the computations predicted well the experimental V s 
values at 0. 1 < Z < 0.2 and S = 7.7 locations. 

Existing experiments 9 measured the velocities at locations 
within very small distance from the wall surface at S = 
7.7, 9.9, and 12.2 stations. The computed V s values at z ~ 
0.01, 0.175 and n = 9600 at locations within N < 0.01 and 
S = 7.7, 9.9 and 12.2 were plotted in Figure 15 and were 
compared with velocity measurements. At z = 0.01, the 
computed V s values agreed very well with the V s 
measurements at S = 7.7 and 9.9. At these S stations, the 
computations predicted high V s at N> 0.0025 and Z = 
0.175 locations. The computed V s was reduced to small 


value at S = 12.2 and the V s profiles showed cross-stream 
inflection point profiles (Figure 1). The measurements 
showed only very small velocity reduction at S = 12.2 
station. The difference could be due to the vortices effect 
on the V s computations. The vortices effects were not 
detected in the experiments. 



Fig 14 Span-wise variations of V s at N = 0.0032 and 
n = 9600 



N 

Fig 15 V s profiles at n = 9600 and S = 7.7, 9.9, and 
12.3 
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CONCLUSION 


SUMMARY 


In an effort to establish a numerical computational 
method to investigate steadiness of turbine blade flow, the 
author used the time-stepping RANS numerical solver in 
the NASA Glenn-HT RANS code to predict the 
steadiness of the boundary layer flows with the effects of 
a concave wall and inlet turbulence. The NASA Glenn- 
HT RANS code, with two-equation turbulence models, 
was established for blade surface heat transfer studies. 
However, this author used a version of One-Equation 
Turbulence Model. This model accounts for the flow 
transition with vorticity at the wall surface. Realistic inlet 
turbulence parameter can also be assigned (see the 
APPENDIX) when the turbulence model is used in CFD 
work. 

As a calibration of the numerical computational methods, 
computations were performed to predict the experimental 
results of the near surface velocity and temperature 
variations in a transitional boundary layer flow along a 
concave wall surface. With specific boundary conditions 
and computational grid, the computations predicted well 
the skin friction measurements at different stream-wise 
locations. The computed tangential and normal velocity 
components, at upstream locations, oscillated periodically 
as numerical time iteration progressed. The velocity 
components also exhibited span-wise periodic variation 
which agreed with the vortices induced upwash and 
downwash at span-wise locations. The velocity 
components at downstream stations did not show the 
coherent structure of vortices effects on the velocity 
components and agreed with the vortices broke down at 
downstream stations. The span-wise temperature 
variation, at downstream station, agreed with wall surface 
temperature patterns from experiments. 

As a test of the present CFD methods to analyze the 
turbine blade flow, the computations were used to 
compute the friction factors and velocity variations at 
near surface locations of a boundary layer flow over a 
concave wall surface with a flat recovery wall surface. 
With specified boundary conditions and computational 
grid, the present computational results of friction factors 
converged to the measurements of friction factors. When 
the convergence was approached, the computed velocities 
at the upstream near surface locations exhibited vortex 
effects on tangential and normal components and the 
components varied periodically in terms of time iteration 
steps. The computations also predicted the measurements 
of the tangential velocity components at locations near the 
junction of the concave wall and flat wall surfaces. 


Wall curvature and inlet turbulence level play important 
roles in flow development over a turbine airfoil. As the 
boundary layer flow develops along a concave wall 
surface, Taylor-Gortler vortices may occur within the 
boundary layer flow and the flow can be steady, 
transitional, unsteady and turbulent. Existing experiments 
showed that combinations of concave wall curvature and 
flow inlet turbulence level induced vortices in the flow. 
Vortex formation and broke up caused transitional and 
unsteady three dimensional flow. This paper presented 
computational fluid dynamic studies of the inlet 
turbulence and wall curvature effects on the flow 
steadiness at near wall surface locations in boundary 
layers. The time stepping RANS numerical solver of the 
NASA Glenn-HT code and a one-equation turbulence 
model were used to perform the computations. 

The present computational methods predicted, with 
reasonable accuracy, the velocity properties at near 
surface locations in a boundary layer flow with the effects 
of concave wall curvature and inlet turbulence. As time 
iteration proceeded, the computed friction factors and 
velocity values converged to their measurements in 
existing experimental studies. Periodic oscillations in the 
computed friction factors and tangential and normal 
velocity components were also found at large time 
iteration steps. The computed velocity components also 
showed the Taylor-Gortler vortices effects on the 
relations among the velocity components. 

The numerical scheme was previously developed for the 
predictions of steady state heat transfer rate of turbine 
blade flows. The present computations used small local 
time step in the numerical scheme and single transport 
equation for turbulence modeling in the computations. 
The scope of the computations was mostly limited to the 
predictions of the near wall surface velocity properties in 
boundary lager flow with concave wall and inlet 
turbulence. However, the computations captured some 
interesting flow phenomena in terms of time stepping 
iteration. It was concluded that the present computational 
methods had the potential to analyze the effects of inlet 
turbulence and concave wall curvature on the steadiness 
of gas turbine blade flow. 
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APPENDIX 


Based on the original One-Equation turbulence 
transport equation 3 , the author wrote 1 the following 
equation, in terms of time, t, and orthogonal coordinates, 
x, y, and z, to model turbulent eddy viscosity, v t . 


^■(p v )+ 

_ d 

p{y + v ) dv 

dx 

a dx _ 

d 

_| 

p{y + v)dv 

dz 

a dz _ 


0 

F 

dy 


d 

dz 


p{v + v)dv 

a dy 


+ pCj\-f t2 )Sv 


p(c wl f w - C bl f t2 /K 2 \v/d) 2 + pf n Au : 


( 1 ) 


The above equation assumes that molecular viscosity v, 
depends only on local temperature and the empirical 

constant, C h , , in the diffusion terms of the original 
version was zero. 


The eddy viscosity, V t , is related to the working 
variable, V , through the following relationships, 

V, = % < 2 > 

f„-pl(p+C,.’) (3) 

and 

X = V/V. (4) 

The production term, S , is given by 

S=S + (v/K 2 d 2 )f v2 

S is the magnitude of the vorticity, d is the distance to the 
wall, and f, 2 = 1 - ^/(l + tf vl ) 


The destruction function, f w , is given by 

/ 

f w ~ § 


r 3 ; \ 1/6 

( 1 + c 6 3 


g ft +C 


w3 J 


where g = r + C w2 (r 6 - r) 
and r = V l{sic 2 d 2 ). 


The transition functions are 

fn =C n g, exp(-C, 2 (ty, 2 / Airjd 2 +g, 2 </ 2 )) 
f t2 =C, 2 exp (~C t4 z 2 ) 

where g / = min(0. 1, An / a> t Ax). 

In the transition functions, (O t is the vorticity at the 

surface and Ax is the grid spacing along the x direction. 
The velocity difference between a field point and the 
surface was All. 

The turbulence model used the following set of empirical 
constants: 


C bl = 0.1355, C b2 = 0.622, a = 2/3, re = 0.41 
C w \ = C hl Ik 2 +(\ + C b2 )l < j,C w2 = 0.3, 
C w3 =2,C vl =7.1,C n =l,C (2 =2,C /3 =l.l 
C =2 

c,,4 z.. 

It was also suggested that x ~ 0.1 be used for free stream 
condition. With the range of x and equations (2), (3), and 
(4), the inlet turbulence level can be specified in terms of 
the molecular viscosity. With x ~ 0.1 and Cvl =7.1 in 
equation (3), eddy viscosity much smaller than the 
molecular viscosity was used for free stream in the model. 

The present author implemented 1 the above turbulence 
model in the NASA Glenn-HT code to model the 
turbulence terms in the 3 -dimensional compressible 
Navier-Stokes equations. 
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